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Abstract 

I reconsider the problem of a raindrop falling through mist, collecting mass, and generalize it 
to allow an arbitrary power-law form for the accretion rate. I show that the coupled differential 
equations can be solved by the simple trick of temporarily eliminating time (t) in favor of the 
raindrop's mass (m) as the independent variable. 
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A perennial homework exercise in differential-equation-based courses in Newtonian me- 
chanics is the problem of a raindrop falling through mist, collecting mass.-'^i'^i^ If the rate 
of accretion is assumed to depend only on the raindrop's current mass (or radius) and not 
on its velocity, then the solution is fairly straightforward. When, by contrast, the rate of 
accretion is taken to depend on both the mass and the velocity, one is faced with a pair 
of coupled differential equations, and the trick for disentangling them can be surprisingly 
difficult to find — not only for the student, but also for the instructor who has forgotten the 
method after some years' absence and must rediscover it (as I can personally testify from 
recent experience) .- 

Here I would like to show that a very general version of the raindrop problem [see ([S]) /(jS]) 
below] can be solved by using a versatile technique that ought to have a place in all students' 
(and instructors') mathematical arsenals: namely, eliminating reference to the old indepen- 
dent variable (here the time t) and temporarily taking one of the old dependent variables as 
the new independent variable. 

Students may remember this trick from the analysis of one-dimensional motion with a 
force that depends on position (x) and velocity (f ) but not explicitly on time: 

m ^ = Fix, v) . (1) 

By using the chain rule 

dv dv dx dv , . 

dt dx dt dx 

we can temporarily eliminate t and instead take x as the new independent variable: this 
yields the first-order differential equation 

dv ^, , 
mv — = r{x,vj (3) 
dx 

for the unknown function v{x). In some cases this equation can be solved explicitly.™ Once 
one has in hand the function v{x), one can reinstate time and solve (at least in principle) 
the first-order separable differential equation 

t = • (4) 

Let us now consider the raindrop problem, which involves a pair of coupled differential 
equations for two unknown functions: the raindrop's mass m{t) and its velocity v{t). The 
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first equation is the Newtonian equation of motion for the raindrop, 

dv dm . . 

m— + V —— = mq , (5) 
dt dt ^ ^ 

which is obtained by the standard procedure of looking at the same collection of water 
particles (the "system") at two nearby times, t and t + At, and writing that the rate of 
change of the system's total momentum equals the total external force on the system. The 
second equation states the hypothesized law of accretion for the raindrop: here I shall 
consider the very general form 

^ = Xm^v^ (6) 
dt ^ ' 

where A > is a constant and a and (3 are (almost) arbitrary exponents. This form includes 
the two most commonly studied cases — namely the easy case (a,/?) = (|,0) [accretion 
proportional to the surface area of a spherical raindrop, with resulting acceleration g/A\ and 
the hard case (a,/?) = (|, 1) [accretion proportional to the volume swept out, with resulting 
acceleration g/1] — but is much more general.- I will show that all these problems can be 
solved by a unified technique. 

First, a few preliminary remarks. Since ([5]) / ([6]) is a pair of first-order differential equations 
for two unknown functions, the general solution will contain two constants of integration. 
Since this system is time-translation-invariant, one of the constants of integration simply 
sets the origin of time. The other constant of integration fixes the relation between the initial 
mass and the initial velocity: that is, it fixes the mass at the moment when the velocity 
has a specified value (or vice versa). The simplest solution arises by demanding that m = 
when f = 0; but we will make some partial progress toward finding the general solution as 
well. 

As mentioned above, all the cases with /5 = are easy: the accretion equation ([6]) can 
be solved immediately for m{t) [it is separable], and the Newtonian equation ([5]) can then 
be solved for v{t) [it is linear first-order with nonconstant coefficients]. The trouble arises 
when /3 7^ 0, as now the equations ([5]) and (El) are coupled. 

To decouple them, we employ the technique mentioned earlier: use the chain rule 

dv dv dm 
dt dm dt ' 

forget temporarily about time, and instead consider the velocity to be a function of the 
mass m (i.e. we temporarily use m as the independent variable). Inserting ([7]) into the 



Newtonian equation ([5]) and using the accretion equation ([6]) to eliminate dm/dt (which 
now multiphes both terms on the left-hand side), we obtain 

v^— + = f-m"" . (8) 

dm m A 

Making the change of variables w = v^'^'^, we find- 

dw 1 + p il + P)9 _„ 

\ w = r m , (9) 

dm m A 

which is a first-order linear differential equation with nonconstant coefficients for the function 

w{m). The integrating factor is m^~^^, and after standard manipulations we obtain the 

general solution^ 



m + 
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_(2 + /?-a)A"" ' mi+/5_ " 
It is most convenient to express the constant of integration C in terms of the raindrop's 

mass mo at the moment its velocity is zero: we have^ 

V = Km^+f [1 - (mo/m)2+^-"]T^ (11) 

where K = [{1 + f3)g/{2 + f3 - a)A]i/(i+^). 

The simplest case is mo = 0. Then we have^^ 

1+/3 



m 



K'vT^ (12) 



where K' is a constant that we need not write out. We now insert f|T2l) into the Newtonian 
equation ([5]). Since each of the terms in this equation is linear in m, the constant K' drops 
out, and we get 

dv 1-a 

So we are done: we have proved that the raindrop falls with constant acceleration ^_^~"^ (yf. 
[When (a, (3) = (|, 0) or (|, 1), this gives 5f/4 or g/7, respectively.] We can, if we wish, plug 
back in to get m{t) = fs:"t(i+/5)/(i-"). 

When mo 7^ 0, the best approach seems to be to insert (fTTj) into the accretion equation ([6]) 
and integrate: 

J m~ W [1 _ (mo/m)2+^-"]"T+? dm = XK^ J dt . (14) 
If /3 = 0, the integral is easy: we get 

m(t) = mo[l + (l-a)At/mJ""]T^ (15) 



vit) = ]-^gt + -^^^ (i _ [1 + (1 _ a)\t/ml-'^r^\ . (16) 



and 

7,m = , 

2 -a" (2-a)A 

If/? 7^ 0, the substitution z = (mo/m)^+^~° allows the left-hand side of ffT4l) to be expressed 
as an incomplete beta function: 

^ mo^5((mo/m)2+^-";a,6) (17) 



2 + p-a 

with a = —(1 — a;)/[(l + /3)(2 + /3 — a)] and b = 1/(1 + /3) [this can alternately be written 
as a hypergeometric function 2^1]. But it seems difficult to make further analytic progress. 
One can in any case see from ( |T^ that the large-time behavior is 

m(t) = K"tT=^ 1^1 + J^akt-'' r=s j (ig) 

and hence 

1-a ( Jl+/3)(2+/?-a) \ 

= 2 + /3-c. ^M^ + E^fe^ 1 > (19) 



fc=i 



from which the coefficients {a^} and {hk} can be determined by substitution into (I5])/(I6]). 
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As one recent article put it, "undergraduate mechanics students are sometimes [my emphasis] 
able to solve the nonlinear dynamical equations of motion to find the deceptively simple ac- 
celeration g/7" [B.F. Edwards, J.W. Wilder and E.E. Scime, "Dynamics of falling raindrops", 
Eur. J. Phys. 22, 113-118 (2001)]. A physics student, answering another student's query on 
an on-line forum, was blunter: "This is a very old problem. Unfortunately, I remember the 
answer, g/7, but I don't remember how you get it. It has an unusual solution. There [is] a 
special substitution that you need to make for the mass, otherwise the problem is insoluble." 
|http : //www . physicsf orums . com/showthread . php?t=198859 

For instance, if F{x,v) = g{x)h{v), then ([3]) is separable [in particular, when F{x,v) = g{x) 
this gives the usual conservation-of-energy equation]. 

If F{x,v) = g{x)v + h{x)v'^, then ([3]) is a first-order linear equation with nonconstant coeffi- 
cients for the function v{x); and more generally, if F{x,v) = g{x)v'^~^ + h{x)v'^, then ([3]) is a 
first-order linear equation for the function v{x)'^. 

Likewise, if F{x, v) = v/[g{v) + h{v)x], then can be turned upside-down to get a first-order 
linear equation with nonconstant coefficients for the inverse function x{v); and more generally, 
if F[x,v) = v/[g{v)x^~^ + h{v)x], then 1^ gives a first-order linear equation for the function 
x[v)^ . 

At least two other cases of ^ correspond to physically realizable (albeit highly artificial) sit- 
uations: namely, {a,/3) = (0,1) and (1,1) arise when the raindrop is constrained (e.g. by a 
massless container) to be a cylinder of fixed base and growing height (resp. fixed height and 
growing base). Note that the cylinder's base can here be of arbitrary shape; it need not be 
circular. 

The case P = —1 needs to be treated separately, yielding 

V = Cm~^ exp 

Of course, all cases /5 < are unphysical. 

The case a = 2 + [5 needs to be treated separately, yielding 

ni + f3)g log(m/mo) \TT? 
^ A mi+/3 ) 

Of course, all cases a > 1 are probably unphysical. 
This assumes j3 > —1. 
This assumes (3 > —1 and a < 1. 



(1 - a)\ 



m 



l-a 



